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ABSTRACT 

The hnearly polarized milhmeter and sub-milhmeter emission in Sagittarius 
A* is produced within 10 Schwarzschild radii of the supermassive black hole at the 
Galactic Center and may originate from a hot magnetized accretion disk, where 
electrons are heated efficiently by turbulent plasma waves. In such a scenario, 
the flux density and polarization are very sensitive to the electron heating rate 
and the inclination angle of disk, respectively, and the major axis of the sub- 
millimeter intrinsic polarization, which is aligned with the rotation axis of the 
disk, is perpendicular to the major axis of the polarized near-infrared emission. In 
combination with MHD simulations, which study the properties of the magnetic 
field and viscous stresses, the current spectral and polarization measurements 
give tight constraints on the model parameters. Simultaneous observations will 
be able to test the model. 

Subject headings: acceleration of particles — black hole physics — Galaxy: center 
— plasmas — radiation mechanisms: thermal — turbulence 

1. Introduction 

Our understanding of physical processes in Sagittarius A*, the compact radio source 
associated with a supermassive [M ~ 3.7 x lO^M©] black hole at the Galactic Center (Schodel 
et al. 2002; Chez et al. 2005a), has improved dramatically since the detection of linear 
polarization of its millimeter and sub-millimeter emission (Aitken et al. 2000). It is generally 
accepted that Sagittarius A* is powered by accretion of the black hole in stellar winds (Melia 
1992; Narayan 1998; Rockefeller et al. 2004). The linear polarization reveals a synchrotron 
origin of the emission and sets strict constraints on the magnetic field, gas density, and 
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electron temperature (Agol 2000; Melia et al. 2000; Quataert & Gruzinov 2000). Recent 
VLBI imaging shows that the 3.5 mm emission is produced within 10 rs of the black hole, 
where ~ 1.1 x 10^^[M/(3.7 x lO^M©)] cm is the Schwarzschild radius for a non-spinning 
black hole (Shen et al. 2005). The emission may originate from a hot magnetized accretion 
disk, and simultaneous spectrum and polarization measurements can be used to estimate 
the model parameters (Bromley et al. 2001). 

High spatial and spectral resolution X-ray and near-infrared (NIR) observations rou- 
tinely detect flares from the direction of Sagittarius A* (Baganoff et al. 2001; Genzel et al. 
2003). The characteristic variation time scale of a few tens of minutes, comparable to the 
dynamical time at the last stable orbit, indicates that they are also produced within a few 
rs of the black hole. The spectra, polarization, and variability, especially when considering 
results from simultaneous multi-wavelength observations, suggest that the NIR radiation 
is emitted through synchrotron processes and the X-rays are likely due to the synchrotron 
self-Gomptonization (Eckart et al. 2006a, 2006b; Yusef-Zadeh et al. 2006a). The flares can 
be triggered by rapid releases of magnetic fleld energy near the black hole, which then heats 
electrons to a few tens of MeV producing the observed emission (Liu et al. 2004; Bittner et 
al. 2007). The long wavelength radio emission is less variable. The low quiescent-state X- 
ray flux and high centimeter radio flux density uncover a non-thermal origin of the emission 
and/or emission from an unbounded flow at large radii (Liu & Meha 2001; Yuan et al. 2002). 
Gorrelated flare activities in the X-ray and radio bands suggest that the radio emission is 
associated with an outflow (Zhao et al. 2004; Yusef-Zadeh et al. 2006b). 

Recent SMA and BIMA observations conflrm the linear polarization of the millimeter 
and sub-miUimeter emission (Bower et al. 2005; Marrone et al. 2006a, 2006b). These 
observations not only show that the flux density and polarization are variable on time scales 
as short as a few hours but also reveal a relatively high mean flux density and hard spectrum 
in the sub-millimeter band as compared with the millimeter spectral bump flrst observed by 
Falcke et al. (1998). The observed position angle of the polarization vector has been used 
to infer an external Faraday rotation measure of ~ —50 rad cm~^ (Macquart et al. 2006; 
Marrone et al. 2007) and a position angle of ~ 170° for the intrinsic polarization, which 
appears to be 90° higher than the position angle of the polarized NIR emission (Eckart et 
al. 2006b; Meyer et al. 2006). 

In this Letter^ we consider a Keplerian accretion flow in a pseudo-Newtonian potential. 
Instead of assuming that the electron and proton temperatures are the same (Melia et al. 
2001), we describe the electron heating by turbulence with a single parameter and treat the 
coohng process more accurately. The ratios of the viscous stress to the magnetic and gas 
pressures can be determined by MHD simulations. These leave three basic model parameters. 
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namely the inclination angle of the disk, accretion rate, electron heating rate. §[2] gives the 
basic equations describing the disk structure. We then discuss the synchrotron emission 
from the disk and its polarization and apply the model to Sagittarius A* in § [Hi In § HI we 
summarize the main predictions, discuss the model limitation, and draw conclusions. 



2. Basic Equations for the Disk Structure 

The basic equations for a Keplerian two-temperature accretion flow in a pseudo-Newtonian 
potential with a magnetic field dominated by the azimuthal component are described by Liu 
et al. (2007). The scale height of the disk is given by 

-rkB{Tp + T,){l + 2f3,y'/' 



H 



GMrUp 



where G, /cb, ""^p, fs = ^GM/c^, c, and Tp denote the gravitational constant, 

Boltzmann constant, proton mass, black hole mass, radius, Schwarzschild radius, speed of 
light, electron and proton temperatures, respectively. Pp = B'^ /STrnksiTe + Tp) is the ratio 
of the magnetic field energy density B'^/Stt to the gas pressure nksiTe + Tp), where n gives 
the gas density. /^From the angular momentum conservation, one obtains the radial velocity 
Vr = —P^PpksiTp + Te)(r — rs)/[mp{GMrY^'^] , where P,y is the ratio of the total stress to 
B'^/Stt (Melia et al. 2001) and we have assumed that there is no angular momentum fiux 
through the disk. So 

GMMml^'^ 

n = - (2^ 

AtiMvWp + T,)]3/2(l + 2(3py/^r{r - rsY ' ^ ' 



The energy conservation equation is given by 



dr 



A 



\kB[T,{a + l + 2Pp)+Tp{2.5 + 2(3p)]+mp 

(3) 

where A is the radiative cooling rate and = (GMrY^"^ / (r — r^) is the Keplerian velocity. 
a = x[3K2,{x) + Ki{x) — 4:K2{x)]/4:K2{x), where x = nieC^/kBTe, Ki refers to the ith. order 
modified Bessel function and rUe represents the electron mass. To obtain the disk structure, 
one also needs to specify the electron heating rate by the turbulent magnetic field (Liu et 
al. 2007; Sharma et al. 2007a). The electron heating time r^c = 3CiH{V(.) /c^^ , where Ci 
is a dimensionless constant, (ve) = 2c{l + x)/x'^K2{x) exp{x) and cs = [ksiTi + Te)(l + 
the mean speed of electrons and speed of fast mode waves, respectively. 

Then we have 

dTe Te ^ Tp-T, A 

dr TacVr TconlVr aukBVr ' 
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where rcoui = 3Trmemp{ve)^/256ne^lnX is the electron-proton energy exchange time via 
Coulomb collisions and In A ~ 15 (Spitzer 1962). 

The cooling is dominated by synchrotron, inverse Comptonization and bremsstrahlung 
processes. In the optically thin regime, the synchrotron and bremsstrahlung cooling rates 
are given, respectively, by 



.n 4:e^ri , o 1. n3a;^ + 12a; + 12 , , 

- (^) - " ^ • 

where (p^), (7b — 1-2, e, and are the mean momentum square of the electrons. Gaunt 
factor, elemental charge unit, and Planck constant, respectively, and all quantities here and 
in what follows are given in cgs units. Most of the thermal synchrotron emission is emitted 
at ue ~ 20z/c = 60eB{x + 1)^ /AirnieCX^ = 8.6 x 10^5(x + l)Va;^Hz (Liu et al. 2006). The 
optical depth at ue can be approximated as te — 3{x + l)HA'^y^c'^ /[8{x + 10)7ik-BTei'^] , 
where the numerical factor is chosen such that the expression is accurate in the relativistic 
regime. When te > 1, the synchrotron cooling is suppressed due to self- absorption: 



^ 8(a: + 10)7rfcBTe4 

3(x + i)cm ■ 



We consider the synchrotron self- Comptonization cooling, then the total cooling rate 



A = Asyn + Aic + Abrem = A,y„(l-87rAsyn^/c5^) ^ + A 



brem 



where Agyn is given by equations ([5]) and ([7]) for te < I and > 1, respectively. 

For given P^, M, Ci, and Te and Tp at an outer boundary, one can solve equations ([3]) 
and (jlj) numerically to obtain the disk structure. The thick lines in the left panel of Figured] 
give profiles of the fiducial model, where (3^ = 0.7, f3p = 0.4, M = 4.0 x lO^^g s-\ Ci = 0.38, 
and Te = Tp = GMmp/bk-Qr at the outer boundary r = lO^r^. For such a low accretion 
rate, cooling is not very efficient and the temperature profiles are determined by the electron 
heating rate. The thin solid and dashed lines give the temperature profiles for Ci = 0.40 and 
0.36, respectively. The other model parameters remain the same. Due to the increase of the 
electron heating rate, the electron temperature is higher for the latter, which has a slightly 
lower proton temperature for energy conservation. The profiles of other quantities do not 
change significantly. We note that the disk structure at small radii is very sensitive to the 
stress at the inner boundary 3rs- Here the stress is chosen such that there is no net angular 
momentum fiux through the disk, which is appropriate for strongly magnetized disks. 
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3. Radiation Spectrum and Polarization 

The formulae for calculating polarized synchrotron emission are given by Melia et al. 
(2001), who treat the flow as a slim disk. Since the millimeter and shorter wavelength 
emissions are mostly produced within 10 r^, where electrons are relativistic, as indicated by 
the profiles of ve and Te, the emission and absorption are dominated by the extraordinary 
emission component and the Faraday conversion dominates the Faraday rotation in the 
millimeter and sub-millimeter bands (Melrose 1997). Then the emission, absorption, and 
dispersion vectors defined by Landi Degl'Innocenti and Landolfi (2004) are in the same 
direction. With the slim disk approximation, the plasma is uniform along a light ray and the 
dispersion effects can be neglected. So the formulae given by Melia et al. (2001) remain valid. 
Considering the vertical structure of the disk and the inhomogeneities caused by turbulence 
may reduce the level of linear polarization and induce circular polarization as implied by 
equations 5.74a and 5.74b of Landi Degl'Innocenti and Landolfi (2004). 

A recent study of the magneto- rotational-instability indicates that (3^ ^ 0.7 when the 
turbulence saturates (Pessah et al. 2006). In principle, (3p can also be determined by 
realistic MHD simulations. Then the observed spectrum and polarization can give tight 
constraints on the disk inclination angle i and Ci. The right panel of Figure [1] compares 
the model predicted spectrum and linear polarization with observations. The thick dashed 
line in the top panel gives the spectrum of the fiducial model with i = 50°, which fits 
the millimeter to NIR spectrum. The spectral index below 100 GHz is ~ 1.7, which is 
determined by the structure of the disk. Emission from different radii peaks at different 
frequencies, above which the source becomes optically thin, and lower frequency emission is 
produced at relatively larger radii. The dotted hne fit the low frequency spectrum with F^, = 
(i//10^'^Hz)°'^ exp[— (i//10^°'^Hz)^/^] Jy. As mentioned in the introduction, this component is 
likely associated with an outflow and its linear polarization is negligible (Bower et al. 2002). 
The thick solid line is the sum of the two, which fits the broadband spectrum. 

Although most of the emission from the disk is produced in the optically thin regime, the 
linear polarization at lower frequencies is suppressed due to the dominance of the unpolarized 
component. The polarization is significant only in the millimeter and shorter wavelength 
range as shown in the middle panel, where the thick line corresponds to the fiducial model, 
and the two thin lines have i = 40° and 60°. The polarization fraction in the sub- millimeter 
range increases dramatically with the increase of i. Our fitting to the observations gives an 
i of ~ 50° ±10°. More accurate polarization observations above 300 GHz may give a better 
measurement of this angle. 

The thin dashed lines in the top panel give the spectra for the two profiles with Ci = 0.36 
and 0.40. The emission spectrum, especially in the NIR band, is very sensitive to the electron 
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heating rate. A good measurement of Sagittarius A*'s quiescent flux in the NIR band can 
lead to an accurate determination of Ci (Ghez et al. 2005b; Clcnct ct al. 2004). The model 
predicts an NIR spectral index of ~ —3.5. Any evidence of a harder quiescent NIR spectrum 
would suggest a non-thermal emission component. 

Observations show that the position angle of the intrinsic polarization is ~ 170° (Mac- 
quart et al. 2006; Marrone et al. 2007). Our model predicts that the emission is polarized 
along the axis of the disk below ~ 1200 GHz. The axis of the disk therefore has a position 
angle of ~ 170°. Future VLBI imaging may be able to test this prediction (Shen et al. 2005; 
Bromley ct al. 2001; Huang ct al. 2007). For small values of the model predicts that 
the position angle of the polarization flips by 90° in the sub-millimeter range due to the 
dominance of the high frequency emission by the blue shifted side. We note that in the NIR 
band, where the emission is also dominated by the blue-shifted side, the corresponding posi- 
tion angle of the polarization is ~ 80°, which is consistent with recent observations (Eckart 
et al. 2006b; Meyer et al. 2006). The bottom panel shows the fitting to the position angle 
of the observed linear polarization. Here an external Faraday rotation has been introduced 
following Macquart ct al. (2006) and Marrone et al. (2007). Such an external Faraday 
rotation may originate from the outflow as suggested by Beckert and Falcke (2002). 



4. Discussion and Conclusions 

We have shown that the disk model can account for the millimeter spectrum and polar- 
ization, the ~ 90° difference in the position angle of the polarized sub-millimeter and NIR 
emissions. Our best-flt model predicts a ~ 8% polarization at 150 GHz, a very soft quiescent 
state NIR spectrum, and a disk inclination angle of 50° with the axis aligning at ~ 170° 
position angle. Future polarimetric spectroscopy and imaging will be able to test the model 
and/or better constrain the parameters. 

Since f3p is currently not well determined by simulations, the fltting to the observed 
data is not unique. For (5^ = 0.7, Pp = 0.1, Ci = 1.08, M = 2 x lO^'^g s^^ and i = 50°, we 
obtain essentially the same emission spectrum and polarization from the disk as the flducial 
model. A better constraint of /3p by simulations is therefore critical to pin down the rest 
of the model parameters with simultaneous spectrum and polarization observations. The 
best-fit orientation of the disk is not sensitive to Pp. However, since most of the emission is 
produced near the last stable orbit 3rs, the stress at this radius and the emission from the 
plunging region may affect the best- fit value of i. A full general relativistic MHD simulation 
will not only address these issues and remove the parameters f3i, and f3p, but also treat 
the vertical structure of the disk and the effects of winds or outfiows properly, which may 
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explain the observed external Faraday rotation and the centimeter emission component (Liu 
& Meha 2001; Yuan et al. 2002; Sharma et al. 2007b). However, to apply MHD simulations 
to Sagittarius A* directly, one needs an advanced treatment of the electron heating by 
turbulence, takes into account the general relativistic light binding effects, and docs the 
transfer of polarized emission in a turbulent flow properly. The plasma dispersion effects 
may induce circular polarization and introduce comparable modification to the spectrum and 
linear polarization as the general relativistic effects (Melrose 1997; Broderick & Blandford 
2004). These investigations may eventually lead to a measurement of the black hole spin. 

Compared with the previous models of a small hot accretion torus (Melia et al. 2001; 
Bromley et al. 2001), we make several significantly improvements. A more realistic treatment 
of electron heating by turbulence replaces the assumption of one temperature flow in the 
original model. This also makes the results almost independent of the outer boundary 
radius and temperature(s) . With the pseudo-Newtonian potential, the inner radius is not a 
free parameter anymore and flxed at the last stable orbit. The zero stress inner boundary 
condition is replaced by the zero angular momentum flux condition. Although other models 
may also explain the observations studied in this paper (e.g. Yuan et al. 2002), the model 
presented here has less parameters and assumptions on the dynamics of the accretion flow 
and/or outflow and therefore likely reveals the dominant physical processes responsible for 
the observed emission. 

This work was funded in part under the auspices of the U.S. Dept. of Energy, and 
supported by its contract W-7405-ENG-36 to Los Alamos National Laboratory. 
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Fig. 1. — Left: Profiles of tlie electron (solid), proton (dashed) temperatures, magnetic field 
(dotted), gas density (long dashed), scale height (dot-dashed), and ue (dot-long-dashed) for 
the fiducial model with (3^ = 0.7, (3p = 0.4, Ci = 0.38, M = 4.0 x 10^^ g s-\ i = 50°, 
and Tf. = Tp = GMrUp/bk^r at the outer boundary. The thin dashed and solid lines show 
the temperature profiles for Ci = 0.36 and 0.40, respectively. The other parameters don't 
change, and profiles of other quantities change very little. Right: Model fit to the broadband 
spectrum and linear polarization. The top panel shows the fitting to the spectrum (Liu et 
al. 2004). The spectra of NIR flares are from recent observations [Gillessen et al. 2006 
(dotted); Hornstein 2007 (solid)]. The dotted line corresponds to an unpolarized component 
with = (z//10^-^Hz)°-3exp[-(z//10^°-^Hz)i/3j jy^ jhe thick dashed line fits the millimeter 

to NIR spectrum with the fiducial model. The thick solid line is the sum of the two. The 
thin dashed lines are for the two profiles with Ci = 0.36 (with a higher flux) and 0.40. The 
middle panel compares the polarization fractions. The data are from Aitken et al. (2000), 
Macquart et al. (2006), and Marrone et al. (2007). The thick line is for the fiducial model. 
The thin lines have i = 60° (with a higher millimeter and sub-millimeter polarization) and 
40°. Bottom: the observed position angle of the linear polarization for the fiducial model. 
The rotation axis of the disk has a position angle of 168°. The solid and dashed lines assume 
an external Faraday rotation measure of -44 and -56 rad cm~^, respectively (Macquart et al. 
2006; Marrone et al. 2007). According to the model, the position angle of the polarization 
flips by 90° in the sub-millimeter range at the frequency, where the polarization fraction 
approaches to 0. 



